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STUDIES OF POSITIVE ION COMPOSITION IN THE 
EQUATORIAL D REGION IONOSPHERE 

by 

R , A . Goldberg and A,C. Aik in 
Laboratory for Planetary Atmospheres 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 

ABSTRACT 

Two daytime D region positive ion composition measure- 
ments at Thumba , India are presented for solar zenith angles 
of 53.2 and 27 . 8°. Comparison of upleg ram with downleg wake 
data shows a large increase in the concentration of heavy 
ions 48 + , N0 +, H 2 0; 55 + , H 3 0 + * (HgO) ; and M + > 65 + for the 
downleg reduced shock condition. Peak concentrations of 
48 and 55 occur at unit optical depth for Lyman alpha 
radiation . The ion 37 + , H^O^'HgO, is dominant for x = 27.8° 

but not for x = 53.2°, consistent above 80 km with an origin 

+ + 

from the X-ray production of 0 2 . It is shown that NO ions 

| | , -j_ 

can be transferred to heavy hydrates 48 , 55 and M >65 
by a reaction chain starting with 

NO + + X t M = NO + • X + M 

where X can be 0 ^ , Ng , C0 2 , or a combination of all three 
depending on the rate of reaction. This chain, together with 

-f. + 

a similar reaction scheme starting with O 2 and ending in 19 , 
37 + and heavier clusters leads to a consistent explanation 
for the hydrated ions observed in the D region. 


INTRODUCTION 


Studies of the positive ion composition below 85 km have 
been conducted via sounding rockets employing quadrupoic ion 
mass spectrometers housed in cryogenical Ly cooled pumping 
systems for several years (eg., Nareisi and Bailey, L965; 
Nareisi, 1966, 1967, 1969; Nareisi and Roth, 1970). The 


results have shown that the D region below 85 km is dominated by 

»+ 


hydrated water cluster ions having mass to charge ratios of 37 , 

^ r- r- + — ~+ , ^ ~ + 


+ , N0 H ’-H 2 0, in- 


H 3 ° ' Kg 0 ; 55 , H^O • (HgO) g ; 19 , H 3 0 ; and 
stead of N0 + as predicted. 

The principal ionization source forming the undisturbed 

D region is the photoionization of nitric oxide (NO) by solar 

Lyman alpha radiation to produce NO 1 . This is supplemented 

o 

by solar X-rays at 2-8 A and galactic cosmic radiation which 
produce 0 2 + and N 2 + (Nicolet and Aikin, 1960). It is necessary 

4. 4 

to find a mechanism which will convert ions such as NO and 0 2 
into hydrated cluster ions. Several mechanisms have been pro- 
posed, e.g. photoionization of CH^ and its subsequent reaction 
with water (Swider, 1970), or the formation and ionization of 
submicronic ice nuclei (Hunt, 1969); Laboratory measurements 
have shown that a reaction scheme exists for the formation of 


the major ion 37 + from 0 + ; this scheme has been applied to 


explain the origin of 37 + in the D region (Ferguson and 
Fehsenfeld, 1969). Difficulties arise because it is assumed 


+ 


that large concentrations of C> 2 are produced by the photo- 


ionization of 0 o ( 1 A ) as proposed by Hunten and McElroy (1968) . 
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Neglect of COg absorption of the ionizing radiation has 
led to a considerable overestimation of the importance 
of this source of ionization in the D region. Huffman 
et hi., (1971). Application of the 0^ reaction scheme 
to NO leads to the production of heavier clusters (e.g. 

48 , 55 and above) rather than 37 , in disagreement with 
earlier measurements. Furthermore, it has been difficult 
to find the necessary reaction to start such a chain 
(Fehsenfeld and Ferguson, 1969) . No satisfactory ex- 
planation has been proposed for the apparent absence of 

-f 

NO coupled with currently accepted values for NO and the 
presence of heavier water clusters. 

In this work we describe positive ion composition 

results obtained from measurements made at Thumba , India 

on March 19, 1970 for solar zenith angles of 53.2 and 

27.8°. The composition was measured using quadrupole ion 

mass spectrometers housed in titanium getter pumped systems 

similar to that described in Goldberg and Blumle (1970) . 

Separation of the payload from the rocket motors permitted 

measurements to be carried out during both upleg ram and 

downleg wake sampling conditions for each flight. Also 

included in the payloads were experiments to measure 

positive ion and electron concentrations, Lyman alpha 

o 

radiations, and the 2-8 A solar X-ray spectrum. The 
results of these experiments are employed herein but are 
reported in detail elsewhere (Aikin et al., 1971; to be 


referred to as Paper 1 throughout the remainder of this 
text) . 

Evidence is presented for the breakdown of anibi.-m 

4 •! 

heavy cluster ions to 37 and possibly 19 by shock wav; 
effects induced by the rocket motion. A qualitative? 
comparison is made between the observed positive ion 
composition measured under weak shock (downleg wake) 
conditions and a theoretical model calculated from re- 
actions involving the transformation of NO* and 0 ^ into 
cluster ions. These comparisons indicate the difficulties 
involved in presenting a quantative account of the posi- 
tive ion distribution at D region altitudes. 

EXPERIMENTAL DESCRIPTION 
A . R ocket Description and Performance 

The rockets, NASA 14.425 and NASA 14.42 4 were Nike 
Apaches supporting nine inch diameter bulbous payloads. 

In addition to an ion mass spectrometer, each payload 

included instruments to monitor solar Lyman alpha 

o 

radiation, solar X-rays between 2-8 A (14.425 only), 
and electron and ion density as described in Paper 1. 

Two axis magnetometers and solar aspect sensors were 
carried for attitude determination. The additional drag 
induced by the bulbous payload is responsible for the 
low apogees obtained. Figure 1 illustrates the performance 
of each payload in terms of Mach number. Figure 2 illust- 
rates the front end configuration of the payload experiments. 


Table 1 provides pertinent characteristics of flight 
performance and launch. 

After second stage burnout, nose cone separation 
occurred followed by payload separation from the rocket 
motor and spectrometer cap ejection. Payload separation 
has two functions. First, it permits less probability 
of contamination from residual outgassing of the second 
rocket stage. Second, it reduces drag effects on the 
experimental system and causes the payload to remain in 
a fixed attitude position (if coning is neglected) to 
altitudes below 65 km on downleg. This feature enables 
direct upleg ram sampling to be compared with direct 
downleg wake sampling in the region of interest. 

B . Spectrometer Characteristics 

Both payloads contained quadrupole ion mass spectro- 
meters housed in titanium getter pumped systems assisted 
by triode ion pumps as described in Goldberg and Blumle 
(1970) . Modifications of that system include two ion 
pumps instead of one to double the efficiency of inert 
gas pumping; an off-center quadrupole rod assembly as 
illustrated in Figure 2 to permit space for additional 
complementary experiments; the increase of rod length 
from three to four inches for improved mass resolution; 
reduction in aperture size to allow less flow and longer 
pumping lifetimes. The spectrometer characteristics are 
shown in Table 2. 
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DATA REDUCTION OF ION COMPOSITION 


Typical raw spectra from 14.124 are illustrated 
in Figure 3. These data are presented with the currents 
displayed on a logarithmic scale (ordinate) in which 
each volt is representative of approximately a full 
decade of sensitivity. Although not shown, the currents 
from 14.425 are found to be somewhat lower than those 
of 14.424, an effect due in part to the difference in 
solar zenith angle conditions but mainly cuased by a 
variation in the mass resolution .settings of the two 
instruments and the downleg operation of the 14.425 
spectrometer at higher internal pressure. The degree 
of resolution is inversely related to sensitivity, and 
accuracy in setting two rocket instruments at the same 
resolution and maintaining such fixed values throughout 
a flight is affected by temperature changes and other 
paramters which may vary between flights. The higher 
pressure operating conditions for the spectrometer on 
14,425 were induced by cap ejection at lower altitude. 

The differences in sensitivity between the two instruments 
are accounted for by the normalization processes. 

The two horizontal ledges on the right hand side of 
each spectrum are results obtained in the high pass filter 
mode of the instrument, i.e. when the D.C. component of the 
field is set to zero. Two values are obtained in this mode, 
the first (A) containing all ions from 1 ^ M < “, the 
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second (B) from 46, 49 < M < », for 14.424 or 14.425 
respectively. Since the sensitivity of the instrument 
becomes much greater in the high pass filter mode, 
measurement of these two values provides a natural technique 
for inflight calibration of this sensitivity. This is per- 
formed by taking the difference between the two values and 
equating it to the sum of currents measured by the spectral 
mode within the difference range. It has also been found 
that the values in modes A and B are not obtained at equal 
sensitivity. The relative sensitivity between these values 
was obtained by studying an E region metallic layer in 
which the only heavy ion constituent of importance was 
iron ( 56 + ) . Since this could also be observed in the 
spectral mode, it was possible to ascertain that total 
mode A was approximately 1/4 as sensitive as total mode B. 
This relative sensitivity was assumed constant in the 
treatment of data at lower altitudes. 

Applying the above relative sensitivity procedures 
to the high pass filter mode data, the composition was 
reduced to absolute values by using the total density data 
obtained from Gerdien and radio absorption experiments 
discussed in Paper 1. All masses were assigned equal 
weight in the normalization process, since the region 
under study is below altitudes where free molecular theory 
plays an important role. Mass dependent flow effects in 




this transist ion region art* not well understood, hut should 
gradually approach a very weak mobility typo dependence 
at the lower altitudes. 

POSITI VE] IO N COMPOSITION RESUL TS 
A . Shock Wave Effect s 

Figure 4 traces the raw current ratio of heavy 
constituents to total constituents during up and downleg 
of 14.424 as obtained from high pass titter mode data. 

An enhancement of this ratio thereby implies a greater 
contribution of heavy constituents to the* total current. 

Of note is the uniform value of this ratio between up and 
downleg in the metallic region (90-100 km) where the major 
contribution to I is iron (56’ t ). The lack of change 

a 

between up and downleg in this region suggests that the 
heavy constituents are either metallic atomic ions and 
hence unaffected by shock effects, or of another origin 
unperturbed by the relatively weak shocks there. Below 
85 km, designated the water cluster region, an enhancement 
of nearly 10 is observed for the ratio, indicating a much 
greater contribution of. heavy water cluster constituents 
in the wake sampling region than in that of the ram. 

These heavier clusters are known to have low energies 
of dissociation leading to the suggestion that shock 
wave temperature and pressure effects are sufficient to 
induce breakup of these heavier clusters to lighter clusters 
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(Narcisi, 1970; Narcisi and Roth, 1970), The evidence 
here appears to verify this suggestion because under ram 
conditions, sampling occurs much closer to that portion 
of the shock where ambient ions have been subjected to 
maximum pressures and temperatures. A certain percentage 
of such ions should also flow into the wake region, but 
a much higher proportion of particles sampled in the wake 
will have traversed across weaker portions of the shock 
at lower pressures, where conditions for dissociation are 
much less probable > Similar effects are also seen xn t’-e 
data of 14,425. 

As shown in Figure 1, 14.425 traveled considerably 
slower than 14.424 in the water cluster domain. Comparison 
of either upleg or downleg data between the two flights 
indicates a higher proportion of heavy clusters (i.e. T D 

D 

in Figure 5 and 6) on the slower vehicle. The effect is 
not as apparent in the downleg individual constituent data 
of 14.425 because of a higher operating pressure (~ 50 times 
greater than on 14.424, cf . Table 2). This was induced 
by a reduced wake effect and lower pumping speed; the 
latter caused by cap ejection at a lower altitude on upleg. 
Data from previous flights on faster vehicles (Mach 6-7) 
through this region indicate that 19 and 37 are the only 
clusters remaining for sampling at such high velocities, 
even under wake sampling conditions. 
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Finally, l lit * height of unit optical depth lor Lyman 
alpha radiation, as measured aboard each payload with an 


EUV detector, is indicated on Figures 5 and 6. This height 
should provide the maximum concentration of' NO . It is 
apparent, especially in the better quality data of Figure 6 
that 48 and 55 peak near this height, suggesting that 
N0^~ has in fact gone into the produc tion <»l these clusters 
(Ferguson and Fehsenf eld , 1969). This additional informa- 
tion provides consistency with the above conclusions, and 
suggests that wake sampled data at lower Mach numbers may 
be sufficiently reliable to represent the* major qualitative 


features of D region composition and clustering. 

B , Qualitative Features and Evaluation 

The normalized positive ion composition profiles are 
illustrated in Figures 5 and 6 for 14,425 and 14.424, 
respectively. The curve marked T„ represents the sum of 

B 

all mass to charge ratios for M > 65. In this work all 
ions are assumed to be singly ionized for identification 
purposes. T^ has been obtained from high pass filter mode 
data using sensitivity adjustments discussed earlier. It 
may thereby be in error by as much as a factor of two or 
three. The sum of all mass currents 46 < M < 66 for 
14.424 and 49 < M < 66 for 14.425 have been deducted from 
I to obtain T D values, in order to prevent double counting 
of masses within these ranges during the normalization 
process. The T curve shown is considered to be the best 
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possible representation of heavy mass ions i rom this uala 
at the present time. 

The sweep rate on both pay Loads provided a height 
resolution close to 1 km at D region altitudes. Some minor 
consti cuent data points are seen to be separated by much 
larger height intervals; missing points are absent because 
of their proximity to the extreme sensitivity limits of 

4. 

the instruments. The absence of NO between 75- and 81 km 
on the upleg data of 14.425 is apparently real. All data 

A. 

shows a strong depletion of NO below the Lyman alpha unit 
optical depth height. In this region the masses do not add 
up to the normalized profile because of the omission of 
sporadic trace constituents (probably of contaminent 
nature in this region of low sensitivity) . Below 81 km, 
the data is sufficiently marginal to be treated on a 
qualitative rather than quantitative basis. 

On both flights, the data cannot be read to accuracies 
better than + 10% because of the logarithmic presentation 
necessitated to gain the broad sensitivity range desired. 
Other possible errors due to unknown corrections for relative 
mass flow, errors in the Gerdien and radio absorption data, 
and ther factors imply error bars approaching + 20% on the 
major constituents and + 100% on the minor constituents, 
the latter value caused by the much stronger contribution of 
noise and other background effects to the small currents 
measured for these constituents. On both shots, the currents 


measured below 75 km are very small, causing measurements 
to be made near the maximum sensitivity limits of the 
instrumentation. The Large values shown are a result 
o! the normalization and thought to be representative* 
of the constituents present on a quaLitativ* basis only. 
However, for the downieg portion of 1.4 42 1, the data may 
be quantitatively reliable to altitudes approaching 65 km. 

On upleg, both flights indicated pour sensitivity up 
to 78 km although cap ejection occurred well below this 
height. On downieg, below the free molecular flow region 
(below 90 km), wake properties of the flow do not appear 
to reduce the quantities of ions reaching the* spectrometer, 
as can be seen in Figure 3. Molecular non-cluster constituents 
(e.g. 0+ N0 + ) are found to be present in the wake with 

t-i 

similar mass ratios to those observed under ram conditions. 

The downieg enhanced values of the heavy cluster constituents 
(e.g. 48 + , 55 + ) are opposite to those* expected from free 
molecular flow or mobility effects, and suggest these 
enhancements to be a lower limit to those obtained with 
the anticipated small corrections from mass dependent 
considerations. The wake data is thereby important in 
providing data down to lower altitudes than obtainable 
from upleg ram, and in establishing the presence of certain 
constituents not observed in the ram configuration because 


of shock effects. Unfortunately, for 14.425, no downieg 
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data of any consequency was obtained below 75 km because 
of a loss of sensitivity related to a higher resolution 
setting coupled with higher operating pressures within 
the instrument. 

From the previous discussion it is clear that modifica- 
tions due to shock effects occur in the sampling of ion 
spectra with rocket-borne ion mass spectrometers. It is 
however, possible to develop a qualitative picture of the 
species and their distribution, A precise quantitative 
analysis will require ion mass spectrometer sampling at 
subsonic velocities. 

The downleg ion density profiles from rocket 14.424 
show the primary ion below 85 km to be 37 with 55 and 
30 + appearing as major ionic constituents. The hydronium 
ion, 19 + , as well as 48 + are secondary. The distribution 
of 55 + traces that of 30 + . In addition, the peaks of 55 + , 
48 + , and 30 + occur at 77 km, which is also the peak of 
production of NO* 1 " (30 + ) by solar Lyman alpha radiation; 
this indicating that these ions originate from this 
ionization source. Above 83 km there is a sudden decrease 
in the hydrated ion cluster densities coinciding with the 
increase in electron density in the lower E region of the 
ionosphere . 

The downleg data of 14.425 are not as good in quality 
as those of 14.424 because of higher operating pressures 
so that care should be taken in reading as much into the 
14.425 profiles. It should be noted however, that in 
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this data 37 is no longer the dominant constituent; 
instead, 30 predominates. Ions of secondary importance 
are 37 + and M f > 65 + . 

INTERPRETATION AND DISCUSSION 

In Paper 1 it was shown that the primary ionization 
sources were solar X-rays, galactic cosmic rays and 
ionization of nitric oxide. Therefore, to explain the* 
observed ion distribution it is necessary to start with 
the ions formed by ionization of N^ and 0^ as well as NO. 
The ion Ng is rapidly converted to 0^ and NO (to a 
certain extent) by the reactions of Table 3A . This implies 
that mechanisms must be found to produce hydrated ions from 
C>2 + and N0 + . 

A reaction sequence has been proposed that begins with 
0 9 + and leads to 19 + , 37 + and heavier clusters (Fehsenfeld 
and Ferguson, 1969; Ferguson and Fehsenfeld, 1969). This 
scheme is shown in Figure 7 with the reactions and rates 
detailed in Table 3B. The scheme requires the formation of 

4. 

0^ and subsequent two-body reactions of the fori 

°4 + + H 2° ^ °2 +,H 2° + °2 

° 2 + ' H 2 ° + H 2 ° * H 3 0+ + ° 2 + 0H 

-> h 3 o + -oh + o 2 

H 3 0 + + H 9 0 + M •* H 3 0 + -H 2 0 + M 


( 1 ) 

( 2 ) 

(3) 

(4) 


14 


•V* sPX 


together with the remaining reactions shown in Table 3B. 
Hydrates with mass greater than 37 are produced when 37 + 
reacts with HgO in a three body reaction. However, the 
predominant ion which must be transferred to hydrates is 

i 

NO . This can be accomplished in an analogous way to 0^ 

formation by first forming the cluster ions N0 + -0g, N0 + *N 2 
+ 

and NO ‘COg. The rate coefficient for the reaction 

NO + C0 2 + Ng -» N0 + -C0 2 + N 2 (5) 

has been measured to be 2.5 x 10""^ cm^sec”^ at D region 
temperatures, Dunkin et al (1971). In addition, they 

have placed upper limits of 5 x 10~ ° cm sec and 

—32 6 —1 + + 

2 x 10 cm sec for NO *N 2 and NO »0 2 formation with 

He and Ar as third bodies but not 0 2 and Ng . Predictions 

based on polarizability by Heimerl and Puckett (1971) and 

Beyer and Keller (1971) would indicate larger rate coef- 

+ 4 * 

ficients for NO °0 2 and NO *Ng then the upper limit set 

by Dunkin et al. The D region ion density data require 

•"32 6 “"X 4* 

a value of 6 x 10 cm sec for reaction between NO and 

its clustering neutral, where the neutral bodies of the 

reaction are each total neutral particle density. The 

mesospheric C0 o concentration is unknown but an estimate of 

-4 

the number density mixing ratio is 3 x 10 , (Hays and 

-33 6 -1 

Olivero, 1970) . This would give 9 x 10 cm sec for the 
equivalent COg clustering rate, but this is insufficient 
to explain the data of 14.424 and would require either a 
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larger rate coefficient or contributions from molecules 
other than CC> 2 . However one must bear in mind that part 
of the problem lies with the data itself, in obtaining 
absolute number densities of specific ambient ions from 
D region ion mass spectrometers at supersonic speeds. 

The ion N0 + C0 o transfers to N0 + *H o 0 by the reaction 

2 A 


N0 + 'C0 2 + H 2 0 -> N0 + -II 2 0 -I C0 2 


( 6 ) 


There is a continued reaction with water until (see table 3A) 


N0 + -(H 2 0) 3 + H 2 0 ■> H 3 0 4 '-(H 2 0) 2 + UNO, 


(7) 


occurs as shown in Figure 7. The ions 19 + and 37* 1 are not 

+ 

formed in this reaction sequence; rather ions such as 55 , 

73 + , and heavier hydrates occur. Since the principal ion 
observed is 37 + , it might be argued that this ion is formed 
by shock wave fragmentation of the heavier hydrates. However, 
there is a difficulty with the reaction scheme in that the 
three body reactions involved in the chain leave most of the 
ions as 48 + rather than as hydrates of greater mass. 

A way out of the difficulty can be achieved by use of 
a reaction postulated by Burke (1970) . The reaction 


NO + * H o 0 + H -» H q O + + NO 
2 o 


( 8 ) 


should be exothermic although it has not been observed thus 
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far in the laboratory. The effect of including this reaction 
i s shown in Figure 8. Here the 37 + distribution is shown 
as a function of rate coefficient for (8) . The distribution 

4 * 

of 73 is also shown. The concentration of hydrated ions 
below 80 km is a sensitive function of (8) , Above 80 km 
37 is controlled by 0 2 and Ng production as shown in 
Figure 8 for the case of no 0 ^ + and Ng + sources. 

In paper 1 it was shown that the "effective recombina- 
tion coefficient", derived from the measured electron density 
profiles and Meira’s (1971) nitric oxide distribution, was 
consistent with the laboratory measurements of the dissociative 
recombination rate of hydrated ions with electrons. These 
rates of reaction are now combined with the Meira NO distri- 
bution and the reactions in Table 3 using the total density 
as the third body, M, to predict the altitude distribution 
of ion composition in the D region. It is assumed that the 
negative ion to electron ratio is less than unity at all 
altitudes above 70 km so that the primary loss process for 
ions and electrons is dissociative recombination. A zenith 
angle of 27.8° corresponding to that of 14.424 is assumed. 

The distribution of neutral atomic oxygen used is approxi- 
mately that of Shimazaki and Laird (1970) . The neutral 
water content is taken to be 4 ppm. Some of the clustering 
reaction rates have an inverse temperature dependence 
(Bohme et al., 1969) which has not been considered here. 


The results of the above calculation aiv illustrated 
in Figure 9, The predicted 37 f distribution is dominant 
but Less in magnitude than that observed The calculations 
predict an increase in the concentration of heavy hydrate 

4, 

ions such as 73 at lower altitudes, not observed in these 
flights because of sweep range limitations but found on 
subsequent flights. Additional computations show that the 
apparent reduction of 37 to secondary status in the data 
of 14.425 is explained above 80 km by the reduced 0 ^ + and 

4. 

Ng production with larger zenith angle. However this 
would not account for the absence of 37 + below 80 km. If 
the Burke mechanism is operative, hydrates would be expected 
at lower altitudes even for larger zenith angles 

The sharp cutoff of 37 + at higher altitudes is governed 
by the reaction 


0 4 + + 0 ■* 0 3 + + Og (9) 

competing at higher altitudes with reaction (1) . It is 
therefore quite sensitive to the concentration of atomic 
oxygen . 

Finally, the observed water cluster light ion distri- 
button (19 , 37 ) under wake sampling conditions must contain 
an unknown percentage of shock dissociated cluster ions in 
addition to the ambient contributions. This implies the 
measured distributions of 19 and 37 to be an upper limit 
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for the ambient distribution of such ions, although more 
nearly correct than values obtained under ram sampling, 
CONCLUSION S 

D region positive ion composition profiles for two 
solar zenith angles (27.8° and 53.2°) have been measured 
and normalized to absolute values with the aid of simul- 
taneous positive ion and electron density measurements. 
Strong evidence for the modification of the ambient 
distribution by the shock wave under ram sampling condi- 
tions is presented, leading to the conclusion that wake 
data is more reliable for proper sampling. The ram shock 

is found to be responsible for the dissociation of heavy 

+ 

water clusters, e.g. HLO 9 (H o 0) , into lighter clusters 

o z n 

dominated by n - 1. Heavy clusters, i.e, n ^ 2, are 
found in much greater abundance in the wake region where 
shock effects due to temperature and pressure increases 
should be less severe. 

Under the assumption that wake composition data is 

more closely representative of the ambient distribution 

than ram sampled data, certain qualitative features are 

+ 

immediately apparent. The 37 cluster ion is found to 
dominate the distribution for x = 27.8°, but reduces to 
secondary importance for x = 53. 2°, ’ suggesting X-ray 
ionization control, at least above 80 km. Water cluster 
constituents 48 + , 55 + , and M + > 65 are found to peak near 


the height of Lyman u unit optical depth, suggesting 
clustering chemistry involving NO* for the formation of 
these constituents. An apparent depletion of NO* below 
this height further enhances this argument. A severe 
cutoff of water clusters is observed to occur between 80 
and 85 km, coincident with a strong positive gradient 
in electron density. 

The wake ion composition profiles together with 
simultaneous measurements of Lyman alpha radiation and 
X-ray flux are used to evaluate ion clustering schemes 
for the conversion of N0 + and 0 ^ into hydrated cluster 
ions. The ion-pair production function is based on neutral 
nitric oxide concentrations consistent with recently reported 
measurements by Meira (1971) , together with ionization of 
0 o and N . It is shown that the 0 9 ^ clustering scheme 
is effective in producing 37 + (H^O 4 " • HgO) above 80 km. 

Concentrations of N0 + , NO^H^O, • (HgO) 2 anc * heavier 

hydrates require a rate coefficient for the clustering re- 
action 


N0 + + X + M -> N0 + - X + M 


( 10 ) 


which is greater than has been measured in the laboratory 
thus far for X being CO^ , N^ or 0^ . It has also been 
necessary to employ uhe reaction 


N0 + °H 2 0 + H •* H 3 0 + + NO 


( 11 ) 
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which is currently not verified by laboratory measurements, 

Much additional laboratory work is necessary before 
clustering can be accepted as the origin of D region 
hydrated ions,, The i'ates of many reactions must be 
determined before theoretical and experimental results 
can agree quantitatively. It is apparent that the D region 
positive ion composition is extremely sensitive not only 
to the primary ions produced (N0 + , 0 ^ , N^ + ) but also 
to minor neutral constituents such as NO, H^O, CO^ 
and II. 

On the experimental side, explicit measurements of 
heavy hydrates will require mass spectrometers sweeping 
to higher mass number under subsonic conditions. Only 
after these steps have been taken can the true positive 
ion composition of the D region be established quantita- 
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Date 
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Launch Time (LMT) 
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FIGURK CAPTIONS 

Figure 1 - Velocity profile of 11,425 and 14.424 in terms 
of Mach number. 

Figure 2 - Front end payload confi gui'a f. i on during flight. 

The spectrometer cap is ejected prior to data samp Lina;. 
Figure 2 - Typical up and clown leg raw spectra obtained from 
the ion mass spectrometer aboard 14.424. 

Figure 4 - Current ratios indicating the relative contribution 
of heavy ionic constituents to the total current sampled. 
Figure 5 - The absolute ion composition distribution from 
14.425 data. 

Figure 6 - The absolute ion composition distribution from 
14.424 data. 

Figure 7 - Thu reaction sequence employed for theoretical 
calculations . 

Figure 8 - The altitude distribution of 37* 1 and 73 4 for 

different rates (k) of the reaction NO^'HgO+H H^O^+NO. 

.f. 

The altitude distribution for 37 when no sources of 

4. 4. 

Og and Ng production are included is indicated by 

-9 3 -1 

o o. This computation is made for k = 10 cm sec 

Figure 9 - Predicted ion distribution for X = 27.8°. 
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